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'Image processing mettJbd and system, and medical examination apparatus, for extracting a 
path following a threadlike structure in an image 



FIELD or TIIC nJ\TTJTIO]^J \ 



The invention relates^tci an image processing method of extracting the points 
of a path following a threadlike strucMre represented in an image. The invention more 
5 particularly relatesJx)-aH-iffiage processmg method of extracting the points representing a 
catheter guide-wire in an X-ray fluoroscopic medical image or thin vessels in an angiogram. 
The invention alsb^relates to an ims^^processing system and to a medical examination 
apparatus such as an X-ray apparatus;^awn& a system and means for image processing. 

The invention is applied to medical imaging systems and to the industry of X- 
ray medical examination apparatus. \ \ \ ^ 

■^T^cegQR euTJD or Tim nn^ENTiON - 
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A method to d|termine anJobject comour, referred to as minimal path, between 
two fixed end points in a 2-D image, is disclosed in ma publication "Global Minimum for 



Active Contour Models : A niiriima^P^ Approac^' by\Laurent D, COHEN and Ron 
KIMMEL, in International Joumal ofCon[T^©E:^sion^24(l), 57-78 (1997). This method 
proposes a technique of boundary detection of objects for shape modeling in 2-D images. 
This method particularly aims at solvmgx^^oimdary detectii^ problem by mapping it into a 
global minimum problem and by determiriing^k'patlTofcminimal length from the solution of 



15 



20 that global minimum problem. The method guarantees ffiat^gl^bal minimum of energy is 
found by minimizing curves between twdend points. This methodVomprises steps of 



manually selecting a start point and an encT 



1 



It in an objfect contouivregion of a gradient 
image ; of propagating a front, in the totality ofs^d^^Si^t^im^^^tarting at the start point, 
in such marmer that this front propagates at lower cost in regions of highygradient values until 



the end point is reached : this determines a cost^fiap-^^^Kch is a totally 

a single minimum ; and a step of back-propagating from the end point t^^w^ds the start point 



by the steepest gradient descent in said totally conv^ 
between the start and end points. 



^ex surface having 
^ards the start point 
surface : this provides a minimal path 
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This publication incli^des by reference a front propagation technique disclosed 
in a publication entitled "A fast marching level set method for monotonically advancing 
fronts" by J. A. SETHIAN in Proc. Nat. Acad. Sci., USA, Vol. 93, pp. 1591-1595, February 
1996, Applied Mathematics. Accordin^to said reference, a front, formed in a 2-D grid of 
potential values, is propagated usi ng aJ ^^st Marching Techniqu e" with a determination of 
the front points. The front is a solution om so-called Eikonal Equation. The Fast Marching 
Technique introduces order in the selectioi^of the grid points and sweeps the front ahead in 
one pass on the 2-D image. The Fast Marchmg Technique comprises marching the Front 
outwards by freezing already visited points aenoted Alive, coming from a set of points 
referred to as Narrow Band, and by bringing new ones denoted Far Away into said Narrow 
Band. The Narrow Band grid^points^^ alwa^s^pdated as those having minimal potential 
values in a neighboring stijucture denote4^Mjn-H^ap and the potential of the neighbors are 
further re-adjusted. 

The method known from COHfiN*s publication constructs the convex surface 

of the cost map using said Fast Marchiiig ted&liqueAwhich provides respectively one path of 

/ ( \ \ 

minimal cost joining the start point to each tespec|ive^oint of the front, said front 

propagating until the end point is\;eached. Tnen, me minimal path is provided by back- 
propagating from the end point to thestaj^pcHm by the Vteepest gradient descent in the 
convex siirface. The numerous pat^s-consteiged by propagating the front forwards and 
joining the start point to the different points of the front fo^ forming the convex surface are 
no longer taken into accoimt. Even^^g^Bg ^j jQ™^g theLtar^point to the end point, in the 
operation~of .forwarding the front, is not mb^teepest gradient descent in the back-propaga tion 



operation. 



So, the final path obtained by this known method does^comprise points 



it com^ri 




extracted by tracking. Neither does 
propagation. 

Besides, it is interesting to note thartfae'^Qix its of a path constructed in the 
operation of marching the front forwards ire points whiih have th^smallest possible 
potentials. Starting at the start point, and goin^^'fecwardslfrom one point to the next point 

^kl Action", i. e. a path on 
which the "Sum" or the "Integral" of potentials calculated over pointV^entials is the 
smallest, though strictly continuously growing as aiunction of the nuqib^er of points present 
on said path between the start point and the current poln^n the front/ 



30 must be at the "minimal cost". So, such a path is a pat 
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A first problem in exmcting a threadlike structure is that said threadlike 
structure may be represented in the original image by a number of thin linear segments which 
are not joined in a strictly continuous n^anner, having "holes" between them, and which are to 
be foxmd among a great number of othenthin xmrelated structures, referred to as false alarms. 
A second problem is that said threadlike structure may be very long and sinuous, so that it 
may be far from a straight line and may GV^n present U-turns along its length, and that it may 
be formed by a great mmiber of points. ^ 

On the one hand, a path const^cted using the front propagation technique 

to solve these problems, due to the fact that 
e. a "Sum" of potentials effectuated along 



described in the known publication is not ada^ 
said front propagation ^^ased orTaJr^ction" 

the constructed path. Befe^^use the t hreadl^-e stnJgture is very long, this Sum of p otentials w ill 



soon become very large on a^pgAfolb 
becomes large, the cost becomes hi gh, and, for i 



readlike structure. When said Sum 
izing costs, the known front marching 



technique may generate a path based on thdnearest^ 

" ^--^ ^lirchine ]e 



alarms in order to follow as few 
points as possible. So, the known front mdrching fec^que may generate a path which is far 
from foUov^ng the sinuous andMbngvthrea^like s^ctujre. 

On the other hand, the minii^ path obtMned by the above described minimal 

may^jiot possiblMprovide extracted points strictly 



ch 



path method is a smoothed path, 
fo llowing a lon ged sinuous threadlike structure 



20 -S UMM.\RY OF TIIE INVENTI Q 





The present invention has for an object to pr oA^de an automatic image 
processing method performed between ^t?^/opredetermined end^oints for supplying a path 
strictly following a long sinuous threadlike structure. The invention particularly has for its 
object to provide such a method of processulg^^mi^gind^ imagA)f such a long and sinuous 
ill-represented threadlike structure and of false alarms and^te ^rovdde extracted po ints of the 
threa dlike structure as a continuous linear stnictog^Sfp ^rts ^not^^d path , suitable for 
improving t he visua l ization of said threadlifeg stru cture. 

It is also an object of the inventionlty^^iQvi^^a^^ examination apparatus 
using an image processing system to carry out this method Wi3^o-pracess medical images. 

Such an image processing method js^clSmed in Claim \\And such a medical 
examination apparatus is claimed in Claim 14. 

jyantages of the method are that it i^utomatic, robust ind reliable ; it 



accurately and securely follows the long sinuous threadlik^"stc ucturg jwithout looking for a 

\ 
I 
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' shorter path and without providin4 holes and false alarms ; it is less calculation-time 
consuming than the known front marching technique, it may be applied to construct 3-D 
images from 2-D data and it is implemented with simple means. An advantage of the 
apparatus is that the visualization is improved of thin long structures in a medical image, such 
as a guide-wire in an angiography image for instance, or brain vessels. 
BRIEF DESCRIPTION OF THE DRAWINGS 
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The invention is described >}iereafter in detail with reference to diagrammatic 
figured, wherein : 

FIG.l is a functional block d\agram illustrating the main steps of the path- 
tracking method ; 



FIG.2A to 2C illustrate the fro 
Marching technique in a grid of points ; 




FIG.3 is Y^mvo oftfe'sye 
a number of points alongsttie track determ 
from the Ancestor or Start Poin^ to th^ J^t 

FIGS.4A and 4B illustratejflie call 



track for determining the curvature of the^ath ; 
FIG.5 illustrate m X 

method. 

DESCRIPTION OF PREFERRED^ 




ropagation according to the Filiation Front 

e assigned to each potential as a function of 
the Filiation Front Marching technique 
hild or current point ; 
tion of the turning angle at a point of the 

having a system for carrying out the 

NTS 




The invention first\elates to an image processing method based on a path- 
tracking operation performed betweeh^two fixed end points for supplying a path strictly 



following a long sinuous threadlike 
In particular, the invention relates t 



^iteld on a background in a digital image, 
image processing method of path-tracking the points 
representing a catheter guide-wire in anX-r^^fiuoroscopy medical image. In another 
particular application, the invention may r^SeTo tB^^a^-trac^g of thin brain vessels 

itton also refeWi ' ^ " 
latibsv ~ ' 



which have a threadlike shape. The inveittion also' relltes'tcr^ sj^stem implementing the 
method and to an X-ray medical examinaubnvapparatusjhaving^uch a system and means for 




image processing and image visualization. 

In the first particular case of gu^6e-wire ^ath-tracki^, the medical image may 
be an X-ray static arteriogram image representn^gat least a blood vessel with a guide-wire. In 
cardiology, such an image may be used to present medical data relaied to the blood vessel for 
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* further medical procedures. The medical procedures using catheter deeply rely on the correct 
visibility of the guide-wire which is ^metallic wire introduced in the vessel for guiding the 
catheter. An image processing operation of path-tracking this guide-wire in an arteriogram 
image, that is the detection and locatioiTyof the points belonging to said guide-wire, can serve 
5 several highly interesting purposes. For example, a binary extraction mask may be 

constructed from the path-tracking operation thus increasing visibility, in order to improve 

\ 

the practitioner's ability to determine medical data. After a complete extraction of the guide- 

\ 

wire points, the guide-wire tip can be located ^d an area of interest may be defined aroimd 
this tip. This enables further local processing\fi)r better visibility enhancement of a tool called 
10 stent introduced in the ves gel for its enlargement. In the second particular case, the medical 
image may be an angiogram of the OTam wtrerAthe vessels have been made as visible as 
possible by injection of a^xtnfr ast flu i^^ the patient. The brain angiogram contains very 
thin vessels which may be very difficult to vi^ual^e. The method according to the inventi on 
permits^ttiep ractitioner to bett^>visuaiiz^m d vessels. In both cases, the data may altemately 
1 5 be usedto"con^^ 3-D images. f j| 

The present path^^ti^kingmperartion uses a Front Marching technique referred 
to as Filiation Front Marching techniquejdenoted FFM, which is not based on a "Sum" or an 
"Integral of potentials" to go from one po^nt to another as known from the state of the art, 
but which instead is base d on tenns^of-^^^^jghted S^jms\^f Potentia ls". It is to be noted that 
20 the "Weighted Sums of Potentials" correspond to calculated terms of "Cumulated Costs", 



\ 



which may not be growdng strictly^qontinuously. Tms property of strictly continuous 
grovmig, which exists in the terms of "Suins^^^o^^ the Prior Art technique, no 

longer exists in the terms of "Cumulated Costs" of the pjesent path-tracking operation. So, 



according to the present Filiation F^JnTMarching technique, the result is that the prope rtyjpf 
25 relying on the^o nstruction of a mimmal path no longer exist s\Actually, the path tracking 
operatioiT]^ on said Filiation Front ^^is hfflg tech nique FF^ aims at prov iding a "Bes t 
Pa^^hich is no LboLund,to-h e the "Minimal Path ". The^Best Path" is defined as the path that 
is validated by using the Filiation Front N^^hingTlT^^ hereafter. 

Now, the FFM technique &^^eomprisel a definition of a fimction of 
30 Cumulated Costs associated to each processed point-,3^j^^^s^ased\on the Potentials of the 
processed points and which is especially adaptddlio^Sfch-tracking lo^g sinuous threadlike 
structures. Such a function formed by terms of Cumulated Costs is nq, ^onger strictly growing, 
thus the cost map which may be constructed fr^Hjjrttiese terms of Cumulated Costs is no 
longer exactly convex, and the required path is no longer allowed to be o|5tained by the 
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* simple steepest descent from the Aid point to the Start point. Using this function of 
Cumulated Costs, the Best Path foimd may not be minimal. Moreover, this Filiation Front 
Marching technique FFM is not based on the known Eikonal Equation used in the prior art, 
but, instead, is based on a specific "Distance" for defining a distance relationship establishing 
a direct filiation between consecutive p^oints which belong to a path and which are denoted 
"Father and Child". 

The path-tracking operatic^ using the FFM technique comprises an initializing 
phase of setting end points for the requirec^Best Path, denoted Start and End Points, which 
are defined in an image of Potentials. Using said relation of direct filiation, in a first 
processing phase of the technique, the front is marched forwards, starting from the Start Point 



which is referred to as the "Ancestor", until iireaches the End Point which is referred to as 
the "Last Child", fflien, in a secbt^Jprocessinglphase of the technique, the required Best Path 
is found by trackingS^^^ardsfror^m^^ Child to the Ancestor. 

Referring toTTGrn which is a fimctional block diagram illustrating the image 



1 5 processing method for path-tra^ing-aJhreaalike structure based on the Filiation Front 



Marching FFM technique/said method comprises :^ 

an operatioml of image /data acquisition from an Original Image denoted OI, 
for instance a medical image such^as abm;diogram irieluding a guide-wire or an angiogram 
including thin vessels to be extracted ; said image dataanclude digital intensity levels referred 



20 to as pixel values, and pixel coordm^^m^the Origmal^Image ; 

an operation 2 of construction of a contrast image IP in which each pixel of the 
Original Image is associated to a n^ calculated intensityAlevel which is taken as a Potential; 



25 



such an operation may be carried ou t by^agy te chipquej known to those skilled in the art for 



enhancing troughs in the OrigtnaKImage : for instance, in the cited medical Original Image 




\ 

OI, the objects of interest, either guide^v^^^opacified vessels, are darker than the 
backgroimd ; so, the Original Image OI is firstmV^rted, then the ridgeness is evaluated by a 
filtering operation : by this filtering operation, fo fiexample, th^maximum curvature is 
calculated at each current point as a\differential invariant by passing a Gaussian filter then a 
differential operator, thus providing enhanced ridges in the p^ce^of the objects of interest ; 
the filtered image is fiirther again inverted to^©\adean image referred to as " Image of 



iterest are again darker than 



Potentials" which is the contrast image IP,^where the ^ 
the backgroimd and form thin troughs substantially contrasting with the background ; 
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an initializing phase 3 of setting end points including a Start point A denoted 

Ancestor and an End point B denoted LastiChild between which a path following a given 

object of interest, that is the threadlike object, is to be determined in the contrast image IP ; 

a first processing phase 4 of p^th-tracking using the firont propagation 

technique referred to as Filiation Front Marching FFM technique, to march a Front forwards 

from the Start Point A to the End Point B ; in ^d first processing phase 4, marching the 

Front forwards according to the FFM technique provides First Tracks for the path ; 

\ I 

a second processing phase 5 of patfi-tracking also using the Filiation Front 



t 



t B considered as the last bom Child and 



ig one First Track backwards ; in order to come 
\ 

ficient to follow the filiation from Child to 



Marching technique 
propagating the front U^kwards 
down said First Track, it is : 

Father and Grandfather until the Ancestor whicl^is the Start point A is reached ; it is different 
from going down the steepe(Sb^::^dient descent of thefknown convex cost map ; the required / 
Best Path is obtained in this second procfcs^g phase^S by said backtracking operation. 

the first processing phase 4 



Using the Filiatio^^Mi£Mai^^ FFM, 



of marching the Front forwards has the advantage o^jiot being time consuming because the 
Front is not wide and does ^ot try many solutions. Th^second processing phase 5 has the 
advantage to provide the retired Best Path having no l^les, including no false alarms and 
strictly following the threadlik^^tructure. The description^f this Filiation Front Marching 
technique FFM includes first a desc^^on^of^^ technique, then a detailed 

description of the fimction of te msjofX)umula te^ ^ ^ 

Referring to FIG.2A to FIG.2C, the Front according to the FFM technique is 



constructed in an iterative maimen^FJG.2A to FIG.2C each represents a grid of points in the 
Image of Potentials IP constructed in ^feitpe^ion^^ the^^riginal Image acquired in the 
operation 1 . This Image of Potentials contains the^tw&'Eiad Points set in the initializing phase 
3, including the Start Point A and the End PointBTAn^a dyan tageous specific Distance law to 



is 



be used in this FFM technique is based'^on the known City Block ^Distance. S aid Dis tance law 
permits of propag atinga_ pn a path constructed 15v| E^^Mtecl mique , only along lines or columns 
of the grid of points, witii one grid point int^al be^^^^^^uc^ssive points referred to as 
Father and Child. Also, an advantageous Criterion of costs, which wdll be further described, 
provides a function of Cxmiulated Costs relat^ to a point of the grid palled Child, said 
Criterion taking into accoimt a function of the Stimulated Costs relajfed to the Father of this 
Child and of the Potential at this Child. 
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Referring to FIG.2A, at a^ertain stage of the front propagation denoted instant 
to, the Image of Potentials IP comprises points already visited which are referred to as Frozen 
Points, located in a zone denoted FZ, andlwhich are represented by dots. These Frozen Points 
are not allowed to change state, that is to have the name of their location FZ or their 



5 associated function of Cumulated Costs modified. This is the reason why they are referred to 
as Frozen Points. The already propagated ^^^^ ^ points including the Start Point A 
which is the Ancestor and Frozen points forming paths going to the last Frozen Points foimd 
before the instant tO. This Frozen zone has an outwards boundary referred to as LOFZ. The 
interior of the Frozen zone FZ are points labeled Frozen Points and having only neighbors 
labeled Frozen Points. 

Besides the Frozen Points, ilhe imafge of Potentials IP comprises points which 
are situated in a zone deriSted't^B^-a^jacgSit the Frozen zone, referred to as Narrow Band 
whose points are represented by crosses^ ^Hiis Narrow Band surroimds the Frozen zone 
external boundary and hasHt;selfa^mer external boundary denoted LONB. So, the extemal 
boundary LOFZ of the Frozen zone is a set of points labeled Frozen Points having at least one 
neighbor labeled Narrow feandT?Sh 

The Image of^^otentials IP^^over comprises points referred to as Far 
Points represented by small squaresTtocated in aszone dbnoted FAR, which includes all 
points extemal to the Narrox^^Ban^^M^^e Frozen^ne\^So, the extemal bovmdary LONB of 
the Narrow Band is a set of jjoints labeled Narrow Baljd joints having at least one neighbor 
labeled FAR Point. 

According to this PEM technique, an^rfperatipn of examining a point of the 
Narrow Band, considered a Child, is pei?of 



e Fathenof this point must be a point of 

the Frozen zone which is a neighbor of the Cmm"and-T^ch is such that the function of 

\ 

Cumulated Costs related to said Cnild associated to said Father is minimal with respect of 
other functions of Cumulated Costs wntSb^my be calculated for said Child associated to 



other possible Fathers. For instance, injJog^ej^mpie^ 
examined. This grid part comprises aV this given ir 

Points Pi, P2, P3^P6^hich \ 
Points P9, P4, P5 whichi 
Points P7, Pg whidr^e^iir 



IG.2A, a part of the grid of points is 



in the Frozen zone FZ ; 

\ ^ 

the Narrow Band NB : 



"I 



These points of the grid are examined and, in the example of FIG. 2^, current point P5 is 
found to present two properties : It is located in the Narrow Band and- its function of 
Cumulated Costs, calculated according to the &;herion which will He further described, is 
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* minimal with respect to those of the omer points of the Narrow Band. The current point P5 is 
then selected as the Narrow Band pointthaving the minimal function of Cumulated Costs at 
this instant tO. According to this FFM technique, said current point P5 having such properties 
is decreed a possible Father. ^ 



A 



5 Referring to FIG.2B, at themext instant tl, in the Image of Potential EP, the 

state of this current point P5 is modifie^ by^arching the limit LOFZ of the Frozen zone 
towards the Narrow Band zone so as Id include point P5 in the Frozen zone. Said limit is now 
referred to as LlFZ^^AirthBf^r*th^om^ P5 is a Frozen point having the function of 
Cimiulated Costs^eviously det^;miined ^^d i^referred to as a Father. Now, according to this 
10 FFM technique, it ft^Qoked^OTthe only Child <^f this Father. 

Referring to FIG.2C, at the instant t2 following tl, in the image of Potentials 
IP, the neighbors of the Fati^r^oijit^^^e examined. According to the Distance law of the 
FFM technique, only admcent points on lines LidVolumns may be neighbors and, unlike the 
method cited as prior ar^the diagonal pointe^lre n^ relevant and are thus not considered. So, 
1 5 according to the Distance^lawof the FFM technique^^in the case illustrated by FIG.2C as an 
example, there are four existing^neig^b^rs^ whichWe. Pg, P45 Pe? P2. AH neighbors must be 
examined and possibly updat ed^-^cui state may be modified, i. e. some of them may become 
Children: ^ 

Neighbors Pe, P2 are in F^. They remain Frozen Points. 

\ / \ 

20 Neighbor Pg^ishiFAJfc: At this instant t2, the limit of the NB zone is 

marched towards the FAR zone so as to include^this neighbor Pg in the Narrow Band zone 

according to a new limit L2NB''oftihrN filing the hole created in the 

Narrow Band at the previous inStanttl by marching the FrJzen Zone limit forwards to reach 

the location LIFZ. The Narrow Band'^ref^abLy^dogs^not remain with a hole. At this instant 

25 t2, point Pg in the Narrow Band becomes a Child whose determined. 

Still referring to Fli^^C^'^^^ln^s^tyil^ and the Cumulated Cost 

of this new Child Pg are determinedSl^lieFather of tie Child P^must be the neighbor which 
is already located in FZ and which has thesinaUg^^^^d^ncti^n of Cumulated Costs. 
Only point P5 answers these conditions t^cau§e^|ias been-^s^cted as having the smallest 

30 function of Cumulated Costs according^to the result of the calci^tion of the Criterion as 

explained in relation to FIG.2B. So this^rozen Point P5 may be decreed to be the Father of 
Child Pg. Now, a new function of Cumulated Costs is to be calculated for this Child Pg fi-om 
said further described Criterion, which takes mto accoimt the function of Cumulated Cost of 



the Father P5 and the Potential at the Child. ' 
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Neighbor P4 1^ in NB. It remains in NB. Until instant t3 this point P4 
of the Narrow Band has had a Father and a related function of Cumulated Costs which have 
been determined at a given previous infant. Let this previous function of Cumulated Costs 
be denoted CCi and this previous Fathen be for example P| of the Frozen zone FZ. For 
searching the Father of this Child P4, the points located in the Far zone are not relevant 
points, said points having infinite costs. Pmnt P2 is in a diagonal position with respect to 

and P5 may be possible Fathers. 



Child P45 so P2 is not relevant. Thus, only P^ 
At the following instant t4, it^ 



s investigated whether the point P5 may be a 



better Father for its^eigEboTP^^an the prd\nx)us Father Pi of said Child point P4 of the 
1 0 Narrow Band. To Qmt end, a new funbt^on pf Cumulated Costs denoted CC2 is calculated for 
this Child P4, takinguitQ^ccount that tW^^inl^>P5 may be a Father of P4. For that purpose. 



this new CC2 takes into accoxirrt-fee^function erf Cumulated Costs related to said possible 
Father P5 and the Potential at said Child P4. Thm'jCCl and CC2 are compared. 



When CC2 > CGirit rtfdS&that it Mnot interesting to make the required path 

15 pass through possible Father P5. So point P5 is notpecreed the Father of P4. 

\ \ / \ 

When CC2 < CCi, it is interesting to make the required path pass through 

Father P5, so said point P5 may b^e-de^reed^flx©^^ Fattier of the neighbor P4 instead of Pi 

and CC2 is the new function of Cumi^aieji Cost relate^to the Child P4 associated to the 

Father P5. / 

20 The Criterion giwig the function of Cumulated Costs CCk related to a given 



current point Pk is a function of the^smn of the minijmrai imong the Cumulated Costs of the 
neighbors of said current point Ptand-ttiHPSl^iaiLatsa^ F|>int. According to the present 
Filiation Front Marching techniqiae FFM, since the cost related to the current point is 
calculated from the Minimum "Min" o£a number of Cxmiulated Costs constituted by the 



25 Cumulated Costs of its neighbors located^t*onie^gri^^^rval on the same line or column, the 
Argument "Arg" of said niinimum^JiMffl"''s^^ o^ the point referred to as 

"Father" of the current point Pk wMeh is the last one^to go out of^the Narrow Band. 



The Father of the current poiht^Pk^isflie last point to go out of the Narrow 
Band in order to become a Frozen Point. The^exG^AR Pomt becomes a Narrow Band Point. 
30 This former FAR point had an infinite Co^ It must be attribute^a new function of 

I \ \ 

Cumulated Cost. The point that is the last one to go out of the Narrovv Band gives its function 
of Cumulated Cost to its neighbor FAR Pointed becomes its possime Father. 

The function of Cumulated Cost GQk may preferably^e determined by the 
Criterion described hereafter providing a robust and reli^^l^value; This 'trit is 
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described in two examples of Techniques referred to as CCk First Calculation Technique and 
CCk Second Calculation Technique, which is a refined version of the First Calculation 



Technique. 



The First CalculatiomTechnique to compute CCk takes into account that the 



5 track between the Ancestor and the current point must have on average the lowest Potential 
value. So this First Technique involves V Potential Mean Value instead of the Potential Sum 
or Integral known from the State of the^Art. According to the Filiation Front Marching 



tween tl 



Technique, the length of the path between the Ancestor and the current point is calculated by 
adding 1 each time a Father fathoms a Q^^^- When the Father is located at the above- 
10 described City Block Distance den otedA])fstan ce Lk from the Ancestor, the Child Pk is located 




at a Distance Lk-t-^wSucKlnaj^^^ 
the current points ar^iupdated, an 
Potential Mean Value (1) 

1 j=k 

T ■ T 

where Qi are the Potentials at the current 
The previous^relation 
(CC,,)(I^> 



en as Lk+i. So, the Distances between the Ancestor and 
\ 

e CCk values are further calculated as follows as a 



(1) 



located between the Ancestor and point Pk. 



cc,= 




written as: 



(2) 



Thus, each time a Child is fatiioniSd^meSDistance ismpdated by adding 1 , and the new 
function of Cumulated Cost is calculated from the previous function of Cumulated Cost, the 
previous Potential and the previous Distance. Unlik^the law known from the State of the Art 
where the paths are associated to ^^mallest 2j9k > which penalizes the long threadlike 

structures, according to the invgntion^said^paths^^^^associated to the smallest Qk, which 
favors said long threadlike structures. 



Referring to FIG:3,^hich shows a curve of ^Ck versus the potential Qk, where 
k is the number of points on the path*=s^arting-frpmJ^j\ncestor, it is to be noted that the 
function CCk is weighted by a constanjUwei^rtJming^Aev^ 1/Lk. Unlike the State of the 



Art, where the ^ Qj^ which is used^ weighted by a weight 4 , according to the invention the 

^""^^^ I \ 

function CCk is weighted by the term 1/Lk, ^vffich^deI)ends on length and makes an 

K is the above-described 

the potential Qk has only a 



advantageous correction of the function evaluatiocL^le Distance! 



(2^ 



City Block Distance. However, it is to be noted that in relation 

\ 

small influence on the CCk value which is a cause of lack of sensftivity to local Potential 
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variations. An advantage is that the)fensitivity of the path tracking technique is substantially 
independent of the path length. 

According to the invenVon, in this First Calculation Technique, the function 
CCic is evaluated using first order recurgive filters which supply the required Weighted Sum 
(2) 

A second Calculation T^hnique is described hereafter for solving said 
problem of lack of sensitivity. Now, irfsaiB Second Calculation Technique, the function CCk 
is calculated usii^g'anaverage effe^^ of p^determmed limited temporal spans defined from 
the current point^swhich permits, as thevpath is progressively constructed, to progressively 
"forget" the data relate4.to^^ii^proces^d i^dor to said temporal spans. So, this Technique 

)a^t info account. 



permits to take Local Events in a given p£ 



\ 



For calculating the function €Ck,\tiie track is searched first only locally in the 
best possible direction, t^i^g-llie'lIocabEvents into account. To that end, only the functions 
of Cimiulated Cost related to the point! found in ^given past are taken into accoimt for 
1 5 calculating the function CCR,^hile thevfimctions of Cumulated Cost related to points which 



where found in a longer past are "forgotten" i. e. they^are not taken into accoimt. This is 
obtained using one parameter a^hi^i^ weightVactor of niininiization of the influence of 
points situated farther away th^n^S^^^^n distanci from the current point. This weight 

factor is applied to the previously described function of Cumulated Costs, where CCk-i is the 

V / ^ 

function of Cumulated Cost related^o the Fathered Qk is the Potential at the current point : 

(3) 

e temporaTspanjWd the number of Fathers which 



CC^ =aCC^_,+(\-a)q-^ 



The a parameter is a constant i 



is taken into account. The span may-.be approximated by the gelation 
l/(l-a) where 0< a < L 




25 If a = 0,9, the number of Fathers which is taken into account to calculate CCk is about 1 0 



If a = 0,5, the number of Fathers wt 



accouQ t^toValculate CCk is about 20. 



Now, it is not sufficient to take Local Events into account, because the 
Filiation Front Marching technique iMy-be.4?£^iv^W^d may lose the good track. So Global 
Events must also be taken into account. To that e^Ffee:^ihationVront Marching technique 
takes into accoimt the curvature along the ^^fSfTra^4n order to avoid zigzags. It is known to 
those skilled in the art to calculate the Turning Angle at a current popit, and to derive the 
curvature from the Turning Angle value. TheTuming Angle is defined as the angle between 



the tangent to the track at the current point and a re|erence axis. T^en a term based on the 

\ 
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'curvature value is taken into accounlfor calculating function CCk in order to penalize track 
trajectories that have too many pointaassociated to important curvature values. 

The Filiation Front Marching technique according to the invention provides a 
front propagation which is made discrete by using the City Block Distance law which permits 
horizontal or vertical increments of value^l , but it does not provide curvature means of 
calculation means. So, according to the invention, a first-order recursive filter is used to 
calculate a reliable curvature along the First Track. The filiation law, based on the City Block 



ion 



Distance, gives the distances DX^, DYk along the lines and the colimms, i. e. the x-axis and 
y-axis, between Fadier Pk-i and^Ghild Bk as :\ 

DXk\[-1,0, 1] y\ \ (5A) 
DYk e [^OrH^ \ \ (5B) 
Diagonal movements are not allowed. So a Tiirnmg Angle is obtained by a calculation based 
on the value of the ratio DXTI^rA^st Turning. 



^\ngle is calculated relative to the short 
relative to the long past. As an example, the 



past including 10 Fathers is^a^^short past (S^^/the past including 20 Fathers is a long past 
(LP). So, the distances are given IS 



past and a second Tumiri^ Angle is calculated relat 

l^rtJieifelations . \ » 
DXsP(K) = /3 DXk + 0^)\dXsp(k-i) \ 
DYsP(K) = /3 DYkMI^^^^ 
Now, taking into account the longoast (LP) to calculate the distances : 
DXlp(k)= rDXK + (K20 
1 D^^= r DYk + (1 - ) u I lp(k-i) 
where 0< /} < y <\, and andy Me-ewKtart^arameter^^Approximately, /? = 0.1 for 10 
Fathers and y = 0.05 for 20 Fathers. From the above relations^of DXk and DYk, 6A, 6B, 7A, 




(6A) 
(6B) 

(7A) 
(7B) 



7B, the Turning Angles are calculated at eaQKE2i2^>5^^^^ taking into accoimt the short past 
SP or the long past LP. These Turning Angles are deni^isd respectively ^sp(K), ^lp(K)- 



\ 



0 sp(K) = atg ( DYsp(K)/DX^(K)) 
^LP(K) = atg ( DYLP(K)/DXLp!lg) 
from which the curvature Kk is calculated at the^ 
Turning Angles relative to short past and longy^ast so as^ 

Kk = I 0 LP(K) - 0 sp(K) I 



active 



ent4)oint by a difference of the two 



(8A) 
(88) 



(9) 

as illustrated by FIGs 4A, 4B showing the possible relative dispositipnsj^of ^sp(k), ^lp(K) 
resulting in different possible values of the dififerenc^v^K- 
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From the curvature Kk, and a new weight W which takes the potential and the 
curvature into account and may be, for inst^ce, W ~ 0,5, the cost Ck is provided by the 
following recursive relation : 

CCk = a CCk-i + (l-a ) [Qk \ W Kk] (10) 
5 According to the above relation, CCk is a funcuon of local measures which provide a good 
sensitivity and of global m^sures^idi pro^^e substantially smooth paths. The calculation 
of CCk may be carried out Bx^ecursiA^fi^re for determining the geometry, i. e. locations, 
and the kinetic, i. e. speed and afeeelgraJiciiC^^^ belonging to a given path. 

Referring to FIG.5, an X-ray medical examination apparatus 150 comprises 
10 means for acquiring digital image datajof^amedical image, and a digital processing system 



120 for processing these data according to the^ocessing method described above. The X-ray 
apparatus comprises an X-ray source 101, a tabiyi02^for receiving a patient to be examined, 
an optical system 103, 104 for suppl^^i^gmag^ to the processing system 120 which has 
at least one output 106 to send image'daS^I^^splay an^or storage means 107. The display 
and storage means may respectively\be the sCTeen 140 and the memory of a workstation 130. 
Said storing means may altemately be extemal storing naeans. This image processing system 



120 may be a suitably programmed computfeikQfthe^orkstation 130, or a special purpose 



IJTs, Memories, Filters, Logic Operators, which 

. ^^-^^ ^ \ 

are arranged to perform the functions of the method steps according to the invention. The 



processor which has circuit means sue 
are arranged to perform the functions < 
workstation 130 may also comprise a keyB 




